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HIGHLIGHTS 


•  The  detailed  ORR  kinetic  mechanism 
on  FeN4-graphene  was  studied. 

•  Both  four-electron  and  two-electron 
ORR  pathways  were  investigated. 

•  The  ORR  favored  the  four-electron 
pathway  on  FeN4-graphene. 

•  The  rate-determining  step  of  entire 
pathway  was  the  reduction  of  02  into 
OOH. 
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Metal-coordinated  nitrogen-doped  carbons  are  highly  active  in  promoting  electrochemical  oxygen 
reduction  reaction  (ORR).  This  study  describes  in  detail  the  ORR  kinetics  on  FeN4-graphene  based  on  a 
density  functional  theory  calculation.  02  molecules  chemisorbed  on  Fe  site  prefer  hydrogenation  into 
OOH  species  rather  than  direct  breakage  of  the  0—0  bond.  The  subsequent  reduction  of  OOH  species  into 
H2O2  has  a  slightly  high  barrier  (1.13  eV).  However,  this  barrier  could  be  bypassed  by  hydrogenation 
dissociation  into  0  and  H2O,  which  displays  a  low  barrier  (0.47  eV).  Further  O  -►  OH  and  OH  ->  H2O 
reactions  are  kinetically  simple.  Throughout  the  entire  ORR,  the  initial  O2  OOH  reaction  determines 
the  total  rate  and  displays  a  reaction  barrier  of  0.62  eV.  This  kinetic  profile  suggests  that  02  molecules  are 
inherently  favorable  for  reduction  into  H2O  on  FeN4-graphene  following  a  four-electron  process. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Oxygen  reduction  reaction  (ORR),  which  is  a  cathodic  half¬ 
reaction,  limits  the  performance  of  fuel  cell  in  fuel  cell  system 
because  of  its  multi-electron  reaction  character  and  low  reaction 
rate  [1,2].  A  high  efficient  fuel  cell  requires  highly  active  ORR  cat¬ 
alysts  to  promote  the  four-electron  reduction  of  O2  into  H2O  [3,4]. 
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Platinum-based  catalysts  are  quite  active  for  ORR,  but  their  prob¬ 
lems  in  durability  and  cost  strongly  limit  their  commercial  appli¬ 
cations  [5,6  .  Considerable  efforts  are  focused  on  developing  new 
affordable  and  stable  catalysts  to  replace  the  platinum-based  cat¬ 
alysts  [6-9],  and  a  variety  of  non-noble  metal  catalysts  are  devel¬ 
oped  [6].  Among  these  catalysts,  transition  metal-coordinated  N- 
doped  carbons  (MNx-C)  are  proven  to  be  quite  catalytically  active 
for  ORR  compared  with  normal  platinum-based  catalysts  and 
exhibit  very  high  stability  [10-23].  Given  these  findings,  MNx-C  are 
potential  materials  to  substitute  platinum-based  catalysts.  In  spite 
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of  significant  progress  on  MNx-C  catalysts,  the  nature  of  the  active 
ORR  catalysis  of  these  catalysts  remains  unclear.  Some  fundamental 
challenges,  such  as  active  catalytic  sites  of  ORR  on  MNx-C  catalysts, 
are  still  on  debate.  Some  authors  supported  that  the  ORR  active 
sites  are  the  M-Nx— Cy  moieties  [24-27].  Other  researchers  pro¬ 
posed  that  transition  metal  itself  does  not  have  a  function  in  the 
ORR,  instead  the  presence  of  transition  metal  center  serves  to 
catalyze  the  formation  of  N/C  catalytically  active  sites  during  the 
pyrolysis  procedures  28-30].  He  et  al.  [31]  showed  that  the  per¬ 
formance  of  CuNx/C  is  significantly  lower  than  that  of  FeN X/C.  Their 
result  indicates  that  the  metal  has  an  important  function  in  the 
electroreduction  of  oxygen.  Bambagioni  et  al.  [32]  obtained  the 
same  conclusion. 

To  further  optimize  the  structure  and  performance  of  MNx-C 
catalysts,  understanding  the  detailed  information  in  the  catalytic 
reaction  thermodynamics  and  kinetics  is  necessary.  Given  that  the 
MNx-C  catalyst  structures  are  slightly  complex,  experimental  in¬ 
sights  into  the  reaction  thermodynamics  and  kinetics  are  very 
difficult.  Theoretical  calculation  is  a  powerful  tool  to  extract 
detailed  surface  reaction  information  at  the  molecular  level. 

Theoretical  calculations  based  on  density  functional  theory 
(DFT)  are  applied  to  MNx-C  systems  for  ORR  [33-35].  Based  on  the 
adsorption  behaviors  of  O2  and  related  intermediates  on  MNx-C 
catalysts,  all  of  the  calculations  resulted  in  MNx-C  systems. 
Graphitic  M-N4  configuration  is  energetically  favorable  and  stable 
in  the  actual  operated  potential  region  [33-35].  Considering  the 
weak  interaction  between  peroxide  and  C0N4— C,  Kattel  et  al.  [33] 
proposed  that  C0-N4  does  not  promote  entire  ORR  and  a  second 
site  for  peroxide  reduction  is  required.  Calle-Vallejo  et  al.  [34] 
showed  that  the  adsorption  behavior  of  transition  metals  is  not 
intrinsic  because  it  can  be  severely  altered  by  changes  in  the  local 
geometry  of  the  active  site,  the  chemical  nature  of  the  nearest 
neighbors,  and  the  oxidation  states.  Kim  et  al.  reported  that  the  Fe 
sites  on  FeN4-carbon  nanotubes  are  stable  adsorption  sites  for  O2 
and  all  intermediate  species  [35],  which  is  similar  to  the  situation  of 
ORR  on  metal-N4  macrocycles  [36-41].  Chen  et  al.  [36]  combined 
experimental  and  theoretical  methods  to  elucidate  the  mecha¬ 
nisms  of  ORR  on  Fe-phthalocyanines  (Pc)  and  CoPc.  They  suggested 
that  the  ORR  follows  the  4e  route  on  FePc  catalysts,  but  the  2e 
route  on  CoPc  was  similar  to  the  results  of  Sun  et  al.  [37]. 

The  aforementioned  studies  mainly  focused  on  thermodynamic 
information,  such  as  the  adsorption  energies  and/or  the  Gibbs  free 
energy  change.  However,  detailed  kinetic  behaviors  of  the  MNx-C 
catalysts  remain  unclear.  These  behaviors  are  important  to  further 
design  and  optimize  catalyst  structure  and  performance.  The  pre¬ 
sent  study  applied  DFT  calculations  to  examine  the  ORR  over  the 
surface  of  a  model  FeN4-graphene  (FeN4-G),  which  emphasized  on 
the  diagnosis  of  the  kinetics  of  the  involved  elementary  reaction 
steps. 

2.  Computational  methods  and  models 

2.1.  Methods 

All  calculations  were  performed  within  the  spin-unrestricted 
DFT  framework  as  implemented  in  DMol3  code  [42,43].  The  Per- 
dew,  Burke,  and  Ernzerhof  parameterization  of  the  generalized 
gradient  approximation  (GGA)  [44]  described  the  electronic  ex¬ 
change  and  correlation  effects  [45].  Double  numerical  plus  polari¬ 
zation  function  basis  sets  were  used  in  the  calculations.  A  smearing 
of  0.002  Ha  to  the  orbital  occupation  was  applied  to  achieve  an 
accurate  electronic  convergence.  A  5  x  5  x  1  k-point  sampling  was 
used.  The  geometry  convergence  tolerance  for  energy  change,  max 
force,  and  max  displacement  were  0.00001  hartree,  0.002  hartree/ 
A,  and  0.005  A,  respectively.  All  atoms  were  allowed  to  relax  during 


geometry  optimization.  The  optimized  C— C  bond  length  of  gra¬ 
phene  was  1.42  A,  which  agrees  well  with  previous  calculations 
[33].  The  adsorption  energy  EadS  is  defined  as  Eads  =  Etot  -  Esur  -  Ex , 
where  Etot,  Esur,  and  Ex  are  the  total  energies  of  the  FeN4-G  with  an 
adsorbed  X  species,  the  FeN4-G  catalysts,  and  the  isolated  X  species, 
respectively.  A  negative  adsorption  energy  indicates  that  the 
adsorbate  molecule  would  be  energetically  favorable  to  be  adduc¬ 
ted  to  the  FeN4-G  catalyst.  The  isolated  X  species  O2,  OOH,  OH,  O, 
H2O2,  and  H2O  were  optimized  in  a  vacuum  box  of 
15  A  x  15  A  x  15  A.  The  optimized  0-0  bond  lengths  of  O2  and  H2O2 
were  1.225  and  1.472  A,  respectively,  which  agree  well  with  pre¬ 
vious  theoretical  0-0  bond  lengths  of  1.22  and  1.48  A,  respectively 
[33].  The  transition  states  were  located  through  the  synchronous 
method  with  conjugated  gradient  refinements.  This  method  in¬ 
volves  linear  synchronous  transit  maximization,  followed  by 
repeated  conjugated  gradient  (CG)  minimizations  and  quadratic 
synchronous  transit  maximizations  and  repeated  CG  minimizations 
until  a  transition  state  is  located. 

2.2.  Models 

This  study  used  a  periodic  (6  x  6)  supercell  of  graphene  with 
FeN4-coordinated  structure,  which  is  more  favorable  energetically 
than  FeN2  and  FeN3  allotropes  [46,47].  Fe,  N,  and  C  atoms  were  all  in 
the  same  plane  in  the  optimized  structure  (Fig.  la).  The  extracted 
structural  parameters  (Table  1 )  were  close  to  previous  calculation 
results  [47],  which  suggests  the  reasonability  of  the  present  model 
and  calculation  condition.  The  ORR  at  solid-liquid  interfaces  was 
complex  because  they  involved  solid  surfaces,  solvent,  and  sub¬ 
strates.  Therefore,  simulating  these  processes  using  a  simple  model 
system  at  the  solid-gaseous  interface  seemed  difficult.  However, 
the  results  of  these  simple  models  contributed  to  a  better  under¬ 
standing  of  the  complex  systems  (heterogeneous  liquid-phase 
catalysis). 

3.  Results  and  discussion 

3.1.  Adsorption  of  O2  on  N-doped  FeN4-G 

All  possible  adsorption  sites  are  considered  for  O2  adsorption. 
The  structural  optimization  shows  that  the  02  molecule  could  only 
be  chemisorbed  on  Fe  with  an  end-on  configuration,  which  agrees 
with  previous  theoretical  results  [35].  This  finding  indicates  that 
the  Fe  site  was  the  active  center  of  the  FeN4-G  catalysts  [19,24-27]. 
In  addition,  the  adsorbed  O2  (Ctyads))  exhibits  two  different  stable 
configurations,  which  incline  its  head  toward  the  formed  five- 
member  ring  (Fig.  lb)  or  six-member  ring  (Fig.  SI).  The  adsorp¬ 
tion  energies  of  the  two  configurations  are  -0.95  and  -0.80  eV, 
respectively;  the  former  configuration  is  considered  in  the 
following  calculations.  After  O2  adsorption,  Fe  atom  protrudes  out 
of  the  graphene  plane  (Fig.  lb)  with  significant  extension  in  Fe-N 
bonds  (Table  1 ),  and  the  0-0  bond  extends  from  1.225  A  to  1.286  A. 
This  result  indicates  a  strong  interaction  between  O2  and  iron  site 
and  reflects  an  effective  activation  of  molecular  O2.  Table  1  also 
presents  the  structural  parameters  and  adsorption  energies  ob¬ 
tained  from  the  adsorption  states  related  to  the  entire  reaction. 

3.2.  Chemical  behavior  of  adsorbed  oxygen  and  formation  of 
adsorbed  OOH 

Under  the  ORR  conditions,  the  two  possible  pathways  for  the 
surface  reaction  of  02(ads)  are  dissociation  into  0(adS)  [Eq.  (1)]  and 
hydrogenation  into  OOH(adS)  species  [Eq.  (2)]. 

02(ads)  -►  20(ads) 


(1) 
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Top  view 


Side  view 


Fig.  1.  Top  and  side  views  of  the  most  stable  configurations  of  FeN4-G  (a)  and  02  molecule  adsorbed  FeN4-G  (b).  Gray,  yellow,  red  and  blue  spheres  represent  C,  Fe,  0  and  N  atoms, 
respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


02(ads)  +  +  e  — ►  OOH(ads)  (2) 

C>2(ads)  dissociation  is  highly  endothermic,  with  a  reaction  en¬ 
ergy  (the  change  of  the  total  energy  between  productions  and  re¬ 
actants)  of  1.02  eV  (Fig.  S2).  The  activation  barrier  of  this  process  is 
approximately  2.53  eV,  which  is  obviously  too  high  to  overcome  at 
the  working  temperature  of  fuel  cell  (approximately  80  °C).  In 
pathway  (2),  02(ads)  can  capture  one  H+  and  e_  to  form  OOH(adS> 
Zhang  et  al.  [48]  suggested  that  O2  can  adsorb  H+  to  form  H+00  in 
acidic  environment  because  the  whole  system  is  neutral.  OOH+ 
could  be  simplified  to  OOH,  and  the  subsequent  adsorbed  H+  could 
be  taken  as  H  by  considering  the  ionization  potentials.  Reaction  (2) 
is  thermodynamically  favorable  (Fig.  2)  because  the  OOH  species 
chemisorbed  on  the  iron  sites  are  relatively  stable,  with  an 
adsorption  energy  of  -1.87  eV.  We  consider  two  different  situa¬ 
tions  dependent  on  the  location  of  H  species  to  investigate  the  ki¬ 
netic  mechanism  of  02(ads)  hydrogenation  into  OOH(adS>  When  H 
species  are  located  at  the  carbon  atom  (Fig.  2a),  capturing  of  H  by 
02(ads)  is  kinetically  difficult,  with  a  barrier  of  1.60  eV.  Alternatively, 
H  species  could  be  stably  co-adsorbed  at  Fe  site  (Fig.  2b),  although  it 
has  a  total  energy  slightly  higher  than  the  former  situation. 
Kinetically,  the  O2  and  H  species  co-adsorbed  on  Fe  site  are  more 
favorable  for  conversion  into  OOH(adS),  with  lower  reaction  barrier 
(0.62  eV). 

3.3.  Chemical  behavior  of  adsorbed  OOH(adS) 

The  reduction  of  OOH(adS)  into  H2O  requires  a  breakage  of  the 
0-0  bond  either  by  one-step  breakage  or  H202-mediated  process. 
The  H2O2  formation  from  OOH  [Eq.  (3)]  displays  reaction  energy 
and  barrier  of  -1.44  and  1.13  eV,  respectively  (Fig.  3a).  The  high 
barrier  indicates  that  the  reduction  of  OOH(adS)  into  H20  is  kineti¬ 
cally  difficult  through  the  H202-mediated  pathway.  This  finding  is 
also  supported  by  the  fact  that  the  0-0  bond  length  in  optimized 
H202(ads)  is  1.960  A,  which  is  markedly  higher  than  that  in  free  H2O2 
molecule  (1.472  A).  The  higher  bond  length  suggests  a  near 
breakage,  which  is  similar  to  the  situations  on  FePc  and  FeP  [36,37]. 
H202(ads)  is  unstable  and  tends  to  dissociate  as  OH(adS)  because  the 
dissociation  process  has  an  extremely  low  barrier  (0.03  eV). 


OOH(ads)  +  H+  +  e-  -  H202(ads)  (3) 

OOH(ads)  — >  0(ads)  4-  OH(ads)  (4) 

OOH(ads)  +  H+  +  e"  -  20H(ads)  (5) 

OOH(ads)  +  H+  +  e-  ->  H20(adS)  +  0(adS)  (6) 


The  one-step  breakage  of  0-0  bond  in  OOH(adS)  is  possible 
through  the  three  processes  expressed  in  Eqs.  (4)-(6).  The  direct 
dissociation  of  OOH(adS)  into  0(adS)  and  OH(adS)  [Eq.  (4)]  exhibits 
reaction  energy  and  barrier  of  0.26  and  1.18  eV,  respectively 
(Fig.  S3).  It  is  endothermic  and  exhibits  a  high  barrier  and  is  thus 
unfavorable.  Alternatively,  the  O-OH  bond  may  break  with  the 
help  of  H,  in  which  the  kinetics  highly  depends  on  the  position  of 
the  introduced  H  species  because  of  the  spatial  asymmetry  of 
OOH(ads>  When  the  introduced  H  species  is  near  the  oxygen  atom 
fixed  at  Fe,  the  dissociation  of  OOH(adS)  with  the  help  of  H  species 
tends  to  produce  OH(adS)  following  Eq.  (5)  but  exhibits  a  high  re¬ 
action  barrier  (1.14  eV;  Fig.  3b).  When  the  introduced  H  species  is 
near  the  free  oxygen  atom  in  OOH(adS>  the  H-assisted  dissociation 
of  OOH(adS)  generates  O  atom  and  H20  [Eq.  (6)].  This  process  is 
kinetically  favorable,  with  a  barrier  of  0.47  eV  (Fig.  3c).  The  formed 
H2O  molecule  drifts  away  from  the  Fe  site  because  of  its  low 
adsorption  energy  (-0.48  eV),  whereas  the  formed  O  atom  tends  to 
remain  adsorbed  at  the  Fe  site  for  further  reduction  because  of  its 
extremely  high  adsorption  energy  (-4.37  eV). 

The  further  reduction  of  0(adS)  into  H2O  should  follow  the  re¬ 
actions  as  shown  in  Eqs.  (7)  and  (8). 

0(ads)  +  H+  +  e^  ->  OH(adS)  (7) 

OH(ads)  +  e  H20(ads)  (8) 

These  two  reactions  are  also  characterized  with  the  H-assisted 
reduction  process,  and  exhibits  reaction  energies  of  -2.24 
and  -1.87  eV,  respectively  (Figs.  S4  and  S5).  This  finding  reveals  a 
highly  exothermic  feature  and  thermodynamic  favorability.  More¬ 
over,  their  reaction  barriers  are  too  small  [0.48  eV  for  reaction  (7) 
and  0.39  eV  for  reaction  (8)],  which  indicates  that  they  proceed  at 
very  high  rates. 

3.4.  Total  reaction  mechanism 

The  relative  energy  and  reaction  barrier  profiles  of  all  reaction 
steps  are  presented  in  Fig.  4,  which  gives  an  overall  picture  of  the 
entire  ORR  on  FeN4-G.  In  this  figure,  the  total  energy  of  the  O2- 
adsorbed  FeN4-G  surface  is  used  as  the  reference  energy  state. 
Based  on  the  energy  of  this  state,  the  relative  energies  of  the  states 
at  the  subsequent  reduction  steps  are  obtained  with  consideration 
of  the  difference  in  atom  composition  48].  The  detailed  data  on  the 
reaction  energy  and  barrier  for  the  different  reaction  steps  are 
summarized  in  Table  2.  According  to  the  calculated  reaction  bar¬ 
riers,  the  reaction  pathway  (as  red  line  in  Fig.  4)  is  the  most 
favorable  pathway  throughout  the  entire  ORR.  In  this  pathway,  the 
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Table  1 

Summary  of  structural  parameters  of  the  FeN4-G  and  related  adsorption  states. 
Note:  £ads.  adsorption  energy  (eV);  h,  the  elevation  of  the  Fe  atoms  above  the  gra¬ 
phene  plane  h  (A);  d,  bond  length  (A);  dN-c(5)  and  dC-c(5),  the  corresponding  bond 
length  in  the  five-member  ring  configuration;  dN-c(6)  and  dc-c( 6),  the  corre¬ 
sponding  bond  length  in  the  six-member  ring  configuration. 


FeN4-G 

02(adS) 

OOH(adS) 

O(ads) 

OH(ads) 

H202(ads) 

H20(ads) 

Fads 

- 

-0.95 

-1.87 

-4.37 

-2.94 

-0.64 

-0.48 

h 

0 

0.372 

0.334 

0.430 

0.346 

0.247 

0.187 

do-Fe 

- 

1.754 

1.784 

1.654 

1.808 

1.796 

2.029 

do-x 

- 

1.286 

1.468 

- 

0.985 

1.960 

0.984 

dw— Fe 

1.902 

1.915 

1.917 

1.935 

1.915 

1.919 

1.908 

dN-c(5) 

1.380 

1.367 

1.365 

1.366 

1.365 

1.364 

1.368 

dN-c(6) 

1.372 

1.372 

1.379 

1.377 

1.377 

1.372 

1.376 

dc- c(5) 

1.430 

1.429 

1.429 

1.430 

1.429 

1.432 

1.433 

dc-c(6) 

1.440 

1.438 

1.438 

1.438 

1.437 

1.438 

1.441 

ORR  proceeds  through  four  major  steps  all  characterized  as  a  hy¬ 
drogenation  reduction  process.  The  02(adS)  chemisorbed  at  active  Fe 
site  is  initially  hydrogenated  as  OOH(adS),  performs  a  hydrogenation 
dissociation  reaction,  and  generates  0(ads)  and  H20(adS).  0(adS)  is 
subsequently  hydrogenated  as  OH(adS)  and  finally  as  H20.  As  indi¬ 
cated  in  Fig.  4,  this  pathway  is  thermodynamically  feasible  and 
exhibits  low  reaction  barriers.  This  result  indicates  that  02  could  be 
readily  reduced  on  FeN4-G  surface  and  produce  H20  following  the 
four-electron  reaction.  The  highest  barrier  (0.62  eV)  appears  at  the 
initial  reduction  step  of  the  02(adS)-to-OOH(adS)  reaction,  which  acts 
as  the  rate-determining  step.  The  relatively  low  value  of  the  barrier 
in  this  rate-determining  step  reflects  that  the  ORR  could  intrinsi¬ 
cally  proceeds  at  a  high  rate,  which  agrees  with  the  experimental 
observations  [10,21  ].  This  event  differs  from  the  situation  of  metal- 
free  N-doped  carbons,  where  the  removal  of  0(adS)  reaction  is  the 
rate-determining  step  [49].  This  phenomenon  indicates  that  the  02 
reduction  could  end  at  the  two-electron  reduction  stage  for  H202 
formation.  Combined  with  the  high  barrier  in  OOH(adS)-to-H202(adS) 
and  the  instability  of  the  H202  adsorbed  on  FeN4-G,  the  optimized 
kinetic  profile  on  FeN4-G  suggests  that  the  four-electron  reduction 
of  02  into  H20  is  inherently  favorable. 
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Fig.  2.  Optimized  structures  for  the  initial  states  (IS),  transition  states  (TS),  and  final 
states  (FS)  of  02  reduction  into  OOH  on  FeN4-G  together  with  corresponding  activation 
barrier  and  relative  reaction  energy,  (a)  02  and  H  species  adsorbed  on  different  sites  of 
FeN4-G.  (b)  Both  02  and  H  species  are  chemisorbed  on  the  Fe  site. 


Fig.  3.  Optimized  structures  for  the  initial  states  (IS),  transition  states  (TS),  and  final 
states  (FS)  of  OOH  reduction  on  FeN4-G  together  with  corresponding  activation  barrier 
and  relative  reaction  energy,  (a)  OOH  reduction  into  H202(adS),  (b)  OOH  reduction  into 
OH(ads),  (c)  OOH  reduction  into  0(ads)  and  H20(ads). 

Electrochemical  reactions  on  catalysts  are  too  complex  to  be 
simulated  by  a  complete  theoretical  description.  The  reaction 
mechanism  may  quantitatively  depend  on  a  number  of  factors, 
such  as  electrode  potential,  electrolyte,  and  adsorbate  coverage. 
Thus,  the  main  features  of  the  ORR  mechanism  may  be  qualitatively 
provided  by  this  simulation.  The  combination  of  experimental 
electrochemical  measurements  (transmission  electron  microscopy, 


Fig.  4.  Calculated  relative  energies  and  activation  barriers  of  full  reaction  pathways  for 
the  ORR  on  FeN4-G.  The  initial  or  final  states  and  transition  states  are  expressed  as  long 
and  short  lines,  respectively.  The  most  favorable  reaction  pathway  is  expressed  as  red 
lines.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 
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Table  2 

The  reaction  energies  (A E)  and  activation  barriers  (£act)  of  different  elemental  re¬ 
action  steps. 


Reaction  step 

A E  (eV) 

£act(eV) 

02(ads)  “ ”  20(ads) 

1.02 

2.53 

02(ads)  +  H(ads)  *  OOH(ads) 

-2.09 

0.62 

OOH(ads)  -*■  O(ads)  +  OH(ads) 

0.26 

1.18 

OOH(ads)  +  H(ads)  - ”  O(ads)  +  H20(ads) 

-2.96 

0.47 

OOH(ads)  +  H(ads)  -*■  20H(ads) 

-2.60 

1.14 

OOH(ads)  +  H(ads)  H202(ads) 

-1.44 

1.13 

H202(ads)  -*■  20H(ads) 

-1.16 

0.03 

O(ads)  +  H(ads)  OH(ads) 

-2.24 

0.48 

OH(ads)  +  H(ads)  H20(ads) 

-1.87 

0.39 

rotating  ring-disk  electrode,  and  so  on)  and  theoretical  modeling  is 
an  efficient  approach  to  the  investigation  of  ORR  mechanisms  and 
identification  of  high-performance  catalysts  for  ORR.  Therefore, 
more  evidence  is  needed  to  support  our  conclusions  in  this  simu¬ 
lation  for  future  experiments. 

4.  Conclusions 

A  DFT  calculation  is  performed  on  the  detailed  kinetic  and 
thermodynamic  behavior  of  the  ORR  on  FeN4-G.  Metal-free  N- 
doped  carbons  are  active  for  ORR  4,8,9],  but  the  O2  molecule  is 
always  chemisorbed  at  Fe  site  for  the  FeN4-G  catalysts.  The  adsor¬ 
bed  O2  is  then  gradually  reduced  into  FI20  through  the  four 
hydrogenation-related  processes:  02(ads)  -►  OOH(ads)  -►  0(ads)  + 
H20(ads)  OH(ads)  +  H20(ads)  2H20(ads).  The  rate-determining 
step  appears  at  the  initial  reduction  step,  02(ads)-to-OOIT(ads)  reac¬ 
tion.  This  finding  suggests  that  the  ORR  efficiency  over  the  catalysts 
is  similar  to  FeN4-G,  which  could  be  increased  by  improving  O2 
adsorption  and  activation.  Throughout  the  entire  ORR  process  on 
FeN4-G,  the  two-electron  pathway  toward  FI2O2  is  unfavorable.  The 
four-electron  reduction  of  O2  into  H2O  could  readily  proceed  by  a 
more  comfortable  kinetic  profile,  bypassing  FI2O2  formation.  This 
finding  shows  the  unique  catalytic  nature  of  FeN4-G  in  the  ORR 
process.  Further  investigations  on  the  ORR  kinetic  behavior  of  MNx- 
C  catalysts  with  different  metals,  N-configurations,  and  carbon 
structures  is  interesting  for  the  design  and  improvement  of  catalyst 
structure  and  ORR  efficiency. 
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